Metric-Dependent Probabilities that Two Qubits are Separable 

Paul B. SlateJ] 

ISBER, University of California, 
Santa Barbara, CA 93106-2150 

(Dated: February 1, 2008) 

Abstract 

In a previous study {Quant. Info. Proc. 1, 397 [2002]), we formulated a conjecture that 
arbitrarily coupled qubits — describable by 4 x 4 density matrices — are separable with an a 
priori probability of w 0.0736881. For this purpose, we employed the normalized volume 
element of the Bures (minimal monotone) metric as a probability distribution over the fifteen- 
dimensional convex set of 4 x 4 density matrices. Here, we provide further /independent (quasi- 
Monte Carlo numerical integration) evidence of a stronger nature (giving an estimate of 0.0736858 
vs. 0.0737012 previously) for this conjecture. Additionally, employing a certain ansatz, we estimate 
the probabilities of separability based on certain other monotone metrics of interest. However, we 
find ourselves, at this point, unable to convincingly conjecture exact simple formulas for these new 
(smaller) probabilities. 
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An arbitrary state of two quantum bits (qubits) is describable by a 4 x 4 density ma- 
trix — an Hermitian, nonnegative definite matrix havingtrace unity. The convex set of 
all such density matrices is 15-dimensional in nature Q . Endowing this set with the 



statistical distinguishability (SD) metric [3], we were able ^n 
(first essentially raised in ^| and later studied further in 



flea 



to address the question 
of what proportion of 



the 15-dimensional convex set (now a Riemannian manifold) is separable (classically corre- 
lated) 0| . The Peres-Horodecki partial transposition criterion [lol.lll1| provides a convenient 
necessary and sufficient condition for testing for separability in the cases of qubit-qubit and 
qubit-qutrit pairs [l^ . 

In |4| we had taken the probability of separability of two arbitrarily paired qubits to be 
the ratio of the SD volume (Vg D ) occupied by the separable states to the SD volume (Vgp n ) 
occcupied by the totality of separable and nonseparable states. We utilized a conjecture 
(combining exact and numerical results) that 



^ 8 



V°+ n = « 5.64794, (1) 

sd 16g0 > v ) 

the full veracity of which has since been formally established, within an impressively broad 
analytical framework jlfli eq. (4.12)]. Then, on the basis of certain extended, advanced quasi- 
Monte Carlo computations (scrambled Halton sequences Q|), used for numerical integration 



in high-dimensional spaces, we had been led to further conjecture in 



a 



that 



TV 6 



VS D = « 0.416186. (2) 

„ SD 2310 v ; 

(In an earlier study jSj, a number of quite surprisingly simple exact results were obtained 
using symbolic integration, for certain specialized [low- dimensional] two-qubit scenarios - 
which had led us to entertain the possibility in |4| of an exact probability of separability in 
the full 15-dimensional setting.) For a scrambled Halton sequence consisting of 65 million 
points distributed over a 15-dimensional hypercube, we had obtained in 4J estimates of 
5.64851 and .416302 of V°+ n and V§ D , respectively. 

Let us note that the interestingly- factorizable denominators in JIJ) and ((21), that is, 1680 = 
2 4 -3-5-7 and 2310 = 2-3-5-7- 11 have been conjectured elsewhere [Uj seq. A064377] (available 
on-line at www.research.att.com/ njas/sequences/Seis.html) to be the two greatest numbers 
(k) for which the sum of the fourth powers of the divisors of k exceeds the fifth power of the 
Euler totient function 4>{k). (4>{k) is the number of positive integers — including 1 — less 
than k and relatively prime to k |l6j.) 



The fundamental conjecture (J2J), of course, directly implied the further conjecture that 
the SD probability of separability is 

Pin = ^ = T^r 2 « 0.0736881. (3) 

V SD 1171 

Now, since the SD metric is identically four times the Bures (minimal monotone) metric jjj, 
we have the immediate consequence that P% U res = ^sd- ( m re g ar d to arbitrarily coupled 
qubit-qutrit pairs, we reached in [12], by similar methods, the conjecture that the corre- 
sponding probability of separability was the ratio (~ .00124706) of 2 20 • 3 3 • 5 • 7 to the 
product of 7r 3 and the seven consecutive prime numbers lying between 19 and 43.) 

The Bures metric plays the role of the minimal monotone metric. The monotone metrics 



comprise an infinite (nondenumerable) class 17], [la, ll9j] , generalizing the (classically unique) 



Fisher information metric 20|. The Bures metric has certainly been the most widely-studied 
member of this class P- Il3 J21I I22I I2 3 . Two other prominent members are the maximal 

25] and the Kubo-Mori (KM) |26l.l27Ll28j| (also termed Bogoliubov-Kubo-Mori and Chentsov 

29]) monotone metrics. 

As to the maximal monotone metric, numerical, together with some analytical evidence, 
strongly indicate that V^+™ is infinite (unbounded) and that P^ ax = 0. (The supporting 
analytical evidence consists in the fact that for the three-dimensional convex set of 2 x 2 
density matrices, parameterized by spherical coordinates [r,8,4>] in the "Bloch ball", the 
volume element of the maximal monotone metric is r 2 sin 9(1 — r 2 ) -3 / 2 , the integral of which 
diverges over the ball. Contrastingly, the volume element of the minimal monotone metric 
is r 2 sin 0(1 — r 2 ) -1 / 2 , the integral over the ball of which is finite, namely n 2 .) 



4, 



In this analysis, we will seek — analogously to our study of the SD/Bures metric in 

- to determine , V^Q and thus their ratio, P^m- (A wiggly line over the acronym for 
a metric will denote here that we have multiplied that metric by 4, in order to facilitate 
comparisons with the results of our previous analysis J^J , which had been presented primarily 
in terms of the SD, rather than the [proportional] Bures metric. The probabilities themselves 

- being ratios — are, of course, invariant under such a scaling, so the "wiggle" is irrelevant 
for them.) In particular, we first find compelling numerical evidence that 

4vr 8 

V K f M = QW SD = - 361 - 468 - ( 4 ) 



Then, using a scrambled Faure sequence j30(] composed of 70 million points, rather than the 



scrambled Halton sequences employed in j^, we obtain an estimate that 

V S KM « 12 - 6822 - ( 5 ) 



This leads us to the further estimates that 



and 



KM 



ps 

KM 



.0350853, (6) 



ps 

1 SD/Bures 



.476147. (7) 



Critical to our analysis will be a certain ansatz that we have previously employed for similar 
purposes in |7] . Contained in the formula for the "Bures volume of the set of mixed quantum 
states" of Sommers and Zyczkowski [3, eq. (3.18)], is the subexpression (following their 
notation) 

Pv T Pfj, 

where p^, p v (//, v = 1, . . . , N) denote the eigenvalues of an N x N density matrix. The term 
(JHJ) can equivalently be rewritten using the "Morozova-Chentsov" function for the Bures 



metric 



13j, eq. (2.18)], 



as 



CBures(p^Pv) = ~ (9) 

Pv + Pfj, 

Qn = (Pv - Pv) 2 C B ures(p^ Pv)/2. (10) 

Our ansatz (working hypothesis) is that the replacement of CBures in the formulas for the 
Bures volume element by the particular Morozova-Chentsov function corresponding to a 
given monotone metric (g) will yield the volume element corresponding to g. We have been 
readily able to validate this for a number of instances in the case of the two-level quantum 
systems [N — 2], using the general formula for the monotone metrics over such systems of 
Petz and Sudar [171 ec L- (3-17)]. 



13l eq. 



We note here that the Morozova-Chentsov function for the Kubo-Mori metric is 
(2.18)] 

logp^-logp^ 

c K m(p^Pv) = • (11) 

Pv- Pn 

To proceed in this study, we first created a MATHEMATICA numerical integration pro- 
gram that succeeded to a high degree of accuracy in reproducing the formula 13, eq. (4.11)], 
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TABLE I: Estimates based on the Bures and Kubo-Mori metrics, using a scrambled Faure sequence 
of 70 million points for numerical integration. In estimating P metr i c , we use the known values of 
^metric S iven © and © — that is > 5.64794 and 361.468 — rather than the estimates of them 
reported in the table 



metric 


ys+n 

metric 


V Ltric 


P s , . =v s /v s+n 

metric " 1 ' 


Bures 


5.64568 


0.416172 


0.0736858 


Kubo-Mori 


360.757 


12.6822 


0.0350853 



for the Hall/Bures normalization constants [31 



for various N. (These constants form 



one of the two factors — along with the volume of the flag manifold [Lj, eqs. (3.22), (3.23)] 
- in determining the total Bures volume.) Then, in the MATHEMATICA program, we 
replaced the Morozova- Client sov function for the Bures metric (jHJ) in the product formula 
(110)1 by the one (fTTj) for the Kubo-Mori function. For the cases, A^ = 3, 4 we found that the 
new numerical results were to several decimal places of accuracy (and in the case N = 2, 
exactly) equal to 2 N<yN ~ 1 ^ 2 times the comparable result for the Bures metric, given by ([12)1 . 
This immediately implies that the KM volumes of mixed states are also 2 N( - N ~ 1 ^ 2 times 
the corresponding Bures volumes (and the same for the KM and SD volumes), since the 
remaining factors involved, that is, the volumes of the flag manifolds are common to both 
the Bures and KM cases (as well as to all the monotone metrics). Thus, we arrive at our 
assertion (jl|). 

We then numerically integrated the KM volume element over a fifteen-dimensional hyper- 
cube using points for evaluation in the hypercube determined by scrambled Faure sequences. 
(As in jj], the fifteen original variables parameterizing the 4x4 density matrices were first 
linearly transformed so as to all lie in the range [0,1].) These ("low discrepancy" [computa- 
tionally intensive]) sequences are designed to give a close-to- uniform coverage of points over 
the hypercube, and thus hopefully yield accurate numerical integration results. 

The results for the two monotone metrics considered so far are presented in Table |U 

8 ^ 
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The tabulated value of PsD/Bures ~ 0.0736858 is closer to the conjectured value of 
0.0736881 than the estimates (based on 65 million points of a scrambled Halton sequence) 
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obtained in whether we scale there by the true value (0) of Vg^ — which gives .0737087 
- or scale by its estimated value there — which gives 0.0737012. (In a fully independent 
analysis based on 70 million points of a scrambled Halton sequence, we obtained estimates 
of VgJ of 362.663, of V£ M of 12.5809 and of P b KM of 0.034805.) 

Associated with the minimal (Bures) monotone metric is the operator monotone func- 
tion, fsuresit) = (1 + t)/2, and with the maximal monotone metric, the operator mono- 
tone function, fmaxif) = 2t/(l + t) jl3, eq. (2.17)]. The average of these two functions, 
that is, faverage(t) = (1 + 6t + t 2 )/(4 + 4t) , is also necessarily operator monotone [18J, 
eq. (20)] and thus yields a monotone metric. Again employing our basic ansatz, we used 
the associated Morozova-Chentsov function (given by the general formula jl?], p. 2667], 
c(x,y) = l/yf(x/y)) 

_ 4(p M + p y ) 
aV£ra9e ~ Pi + Gprfu + PV [ j 

in the same quasi-Monte Carlo computations based on the scrambled Faure sequence 

composed of 70 million points. We obtained estimates of V^^: age ~ 1.00888 • 10 35 , 
^average ~ 2.825 • 10 22 , so it appears that their ratio, Pa VeTage , is quite close to zero, if 
not exactly so (as we reasoned for the case of the maximal monotone metric itself). 

Additionally, in an independent set of analyses (Table HI}) , based on 160 million points of a 
scrambled Halton sequence, we sought to obtain estimates of the probability of separability 
of two arbitrarily coupled qubits based on three other monotone metrics of interest. These 
correspond to the operator monotone functions, 



l/./r , i\2. f u\ j.t/(t-i)/„. f u\ 2(£ 1 



<2 



fwy(t) = -(Vt + l) 2 ; f G Ks(t)=t^/e; f NI (t) = - > (14) 

4 (1 + £)log(£) 



The subscript WY denotes the Wigner-Yanase information metric 



33, sec. 4] 



the 



subscript GKS, the Grosse-Krattenthaler-Slater (or "quasi-Bures" ) metric (which yields the 



common asymptotic minimax and maximin redundancies for universal quantum coding 



sec. IV.B] 



35 



36j), and the subscript NI, the "noninformative" metric (termed the Morozova- 



Chentsov metric in 37]). 

On a concluding note, let us indicate that the seventy million points of the scrambled 
Faure sequence, central to our analyis, were computed in seven blocks of ten million each, on 
seven independent / non-communicating parallel processors (initialized with different random 
matrices). It is not totally clear to us at this stage, whether or not this is likely to be 
inferior to a (more time-demanding) computation on a single processor, or to the use of 



TABLE II: Estimates based on the Wigner-Yanase, Grosse-Krattenthaler-Slater and Noninforma- 
tive monotone metrics, using a scrambled Halton sequence consisting of 160 million points. In 
estimating P metric , as a ratio, we are compelled to use the estimated values of V^^ ric , unlike Table 
I, since their true values are not presently known 



metric 



ys+n 

metric 



V s 



metric 



P s . . = V s iv s+n 

metric " I " 



WY 
GKS 
NI 



446.615 
166.906 
3710.31 



2.1963 
.9938 
12.616 



.0504375 
.0610692 
.0348745 



seven processors all initialized with identical random matrices (cf. |38|). (In any case, it is 
certain that no such question needs to arise for the [fully deterministic] scrambled Halton 
sequences, also used here and in ^|.) This matter is presently under our investigation. 
Nevertheless, in any case, we have achieved here rather impressive, improved convergence 



to our previously conjectured exact value Q of P; 



SD/Bures- 
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